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Abstract 

Within the framework of the Projective Unified Field Theory the distribu- 
tion of a dark matter gas around a central body is calculated. As a result the 
well-known formulas of the Newtonian gravitational interaction are altered. 
This dark matter effect leads to an additional radial force (towards the center) 
in the equation of motion of a test body, being used for the explanation of the 
so-called "Pioneer effect", measured in the solar system, but without a con- 
vincing theoretical basis up to now. Further the relationship of the occurring 
new force to the so-called "fifth force" is discussed. 
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1 Idea of dark matter accretion around a center 

Application of the author's 5-dimensional Projective Unified Field Theory (PUFT), 
being published in a series of papers (Schmutzer 1995a, 1995b, 1999, 2000a, 2001), 
to a closed homogeneous isotropic cosmological model for the whole time scale 
(big start to presence) led to following numerical results (Schmutzer 2000b, 2000c), 
where the index p means present and y year: 



a) 


tp 


= 18 • 10^ y 


(age of the cosmos), 




b) 




- 77.6 

sMpc 


(Hubble factor). 
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c) 




= 3.3 • 10^27 g(,j^-3 


(mass density). 





Empirical astrophysical estimates of the visible mass density of the normal 
(mostly baryonic) matter come to the maximum value ss 10"'^'^ gcm^^. We in- 
terpreted our result with respect to this difference in the mass densities as a hint 
at the existence of a dominating dark matter part with a mass density more than 
two orders of magnitude greater than that of the normal matter. 

Our hypothetical model of the cosmos investigated, roughly corresponds to the 
following picture of the present cosmos: 

Our cosmos consists of a gas of dark matter particles (dm-particles) penetrating all 
matter of the cosmos, particularly also the existing compact objects (stars, nuclei 
of galaxies etc.) which look like buoys in the dark matter sea. Assuming a homo- 
geneous gas of one sort of particles, by some hypothetical arguments we recently 
were led to the following numerical values (Schmutzer 2000b): 

a) nip = 2 • 10~^^ g (rest mass of the dm-particle), 

(2) 

b) Tp = 1.78 K (kinetic temperature of the dm-gas). 
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As a result of the attractive interaction of the dm-particles with the compact 
object considered we expect a statistical distribution of the dm-particles within and 
around the compact object being embedded in this dm-gas. Let us introduce the 
notion "accretion cloud" (Umwolke) for this part of the dm-gas which by accretion 
exceeds the cosmological average distribution of the dm-gas. For simplicity we 
model the compact object by a homogeneous mass sphere. Further for simplicity 
we apply the Boltzmann-Maxwell statistics describing spinless particles without 
degeneracy. 

Let us here mention that in a different context and on a different theoretical 
basis H. Dehnen et al. (Dehnen 1995) used the Bose-Einstein statistics (preferred 
by physical arguments) and the Fermi-Dirac statistics (dropped) in treating the 
dark matter in galaxies in order to find a theoretical concept in understanding the 
problem of the "flat rotation curves" of stars in galaxies. In interesting papers 
detained quantum-statistical calculations were performed. We think that for a 
rough understanding of our approach it is legitimate to simplify the situation by 
referring to the Boltzmann-Maxwell statistics. 



2 Basic equations 

In the following PUFT is specialized to the nonrelativistic and weak field case, but 
we take into account the scalaric effects described by the scalaric field function a. 
Then the gravitational field equation is given by 

A$ = ATTjNfJ' (3) 

(<i>(r,t) local Newtonian gravitational field of the central body, fi mass density of 
the central body, jn Newtonian gravitational constant as a true constant). 
The scalaric field may consist of two parts according to 

a) a^ac + s with b) |s| < |(Tc| , (4) 

where adt) is the global scalaric cosmological field and s{r,t) the local scalaric field 
of the central body. Then the scalaric field equation reads 

Act = As -— . (5) 

In the general case of a test body moving in the external fields $ and s as well as in 
eventually existing external electromagnetic fields E and B (neglecting in all fields 
the back reaction) the equation of motion reads: 

M(^+gr.d^ + ^Sr^ds + v^)=Q(E+^) (6) 

\dt (Tc At J \ C J 

(M inertial mass, Q electric charge, v velocity of the test body). Let us mention that 
the last term on the left hand side of this equation of motion leads to "bremsheat 
production" , particularly also in celestial bodies, investigated in detail in previous 
papers (Schmutzer 2000c, Schmutzer 2001). 

In the following we refer to a test body without electric charge, i.e. Q = 0. 
Because of the same mathematical structure of both the equations (^ and (^, in 
this approximation it is possible to introduce the combined field function 

* = $(^1-^^ +c2lna, + ^$, (7) 

($c(i) cosmological value of the Newtonian potential). Hence we obtain the common 
field equation 

A*^47r7Ar (l-^)/^ (8) 
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and the equation of motion 

dii + g,ad* + ^^^0. (9) 
Let us now use following notation: 

fl = flri+ tJ-dm (10) 

(/i„ mass density of normal matter, ^dm mass density of dark matter). Application 
of this latter formula to the interior (i) and to the exterior (e) of the spherical 
central body leads to the equations 

a) = Mm + Mdmi and b) = Mne + Mdme ■ (H) 

Next we split the mass densities into both the parts coming from the global cosmo- 
logical matter (normal matter and md-matter) and from the local matter (central 
body and accretion cloud): 

a) l^ni — l^ng ^" f-^ns: b) l^dmi — l^dmg ~t~ ^dmi (-^^) 

and 

a) l^ne — l^ng ; b) l-^dme — lading ~t~ f^dme ■ (-^'^) 

Here the indices refer to: 

g — > global (cosmological contribution), 
s — > sphere (central body contribution). 

The roof denotes the dm- contribution induced by the central body, i.e. surplus 
dark matter to be treated statistically. 

For our further calculations it is convenient to introduce the abbreviations 

a) Me = MG + Arfme with b) flG = Mng + Mdmg , i-B- 

(14) 

C) Mi = MG + Mns + fj-dmi ■ 

Let us remind that our cosmological model led us to following numerical values of 
the scalaric cosmological world function: 

a) crc(i « 1000 y) « 3.3, b) crcp = 65.188, i.e. 

(15) 

c) ^ 1 (for the time scale under consideration). 
Then the field equation (^ for the interior and exterior reads: 

a) A*, = 47r7Ar/i,, b) A-^e ^ in-fNjle- (16) 
Further from (|^) results 

^ = <P + c^\nc7c- (17) 

Next, analogously to the above formulas for the mass density of the dm-particles 
we now write down the corresponding formulas for the particle number density ndm 
being related to the mass density as follows (m mass of a dm-particle): 

a) /idm = mUdm , b) ^J,dmg = mUdrng , c) fldm = mfldm ■ (18) 

We find 

(19) 

In ( p^ ) and (|l^) we omitted the idices i and e. 

Since the time dependence of the cosmological quantities (e.g. of adt) ) is 
extremly slow (adiabatic time dependence), compared with the time dependence 
of the cosmogonic processes of the cosmic objects, we omit explicit writing of the 
time t in such cosmological quantities. Particularly constants of spatial integration 
implicitly involve this adiabatic time dependence. 
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3 Statistical dm-particle distribution in the accre- 
tion cloud 



3.1 Distribution function 

According to our accretion concept presented in Section |l| the statistical dm-particle 
number density for the case of spherical symmetry reads (k Boltzniann constant): 



ndr, 



exp 



kT 



1 



(20) 



where n denotes the dm-particle number density far away from the cental body: 
a) fi — n{r = oo) with b) Xoc = Xe{f = oo) = 

The potental function x used is related to by 
a) = X + *c , where b) ^'c = 

Using (pn) we can write the equations (Iq) as 



In <7c 



(21) 
(22) 



a) = 

b) Axe = 
where the quantity 

Mo — MG + Mris 



47r7jv 
47r7jv 

— const 



^0 + mn |exp(-^) - l| 
+ run |exp(-^) - l| 



(referring to a homogeneous sphere) means the mass density of the sphere. 
Let us for the following use the approximation assumption 



mx 



kT 



< 1 



being well fulfilled in our further applications. Hence series expansion of i 
to the linearity 

nmx 

n = 

kT 



(23) 
(24) 

(25) 
leads 

(26) 



between n and x- This simplification means that we are able to perform the further 
calculations analytically. Then the equations (|2^ ) commonly treated (omitting the 
indices i and e) read 

Ax + K^x - 47r77vA , (27) 
where /t means /io for the interior or fiQ for the exterior, and is defined by 



kT 



(28) 



Let us here mention that the author's dissertation had referred to the theory of 
strong electrolytes, where according to the Debye-Milner theory the statistical treat- 
ment of a (negative) ion cloud around a fixed (positive) ion (and vice versa) plays an 
important role. This knowledge induced the author's idea to treat the dark matter 
gas analogously. But in the first moment he was surprised that, though between 
positive and negative electric charges as well as between masses attractive forces 
act, in his basic equation ( p7| ) the second term on the left hand side appeared with 
a positive sign, in contrast to the theory of strong electrolytes. The consequence of 
this fact is that in the Debye-Milner theory a Yukawa-like term with exponential 
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decrease occurs, whereas in this theory an interesting phenomenon of periodicity is 
met. 

In solving the equation ( |27| ) we first go over to spherical polar coordinates: 

^(r'^)+^'r\^An^^fir\ (29) 
dr \ dr J 

3.2 Interior 

In this case the substitution /i — > /ip has to be applied. Then equation ( |2l| ) reads 

A (r^^) + = 4^7^Mr2 . (30) 

dr \ dr J 

The solution of this inhomogenious differential equation is 
B B 

Xi — — cos(/vr) H sin(Kr) + Bq (31) 

r r 

a) Bq = — — = — , where b) Mc = — ^ (32) 

is the mass of the central body with radius tq (-Biand B2 constants of integration). 
Regularity at r = leads to _Bi = , i.e. 

B2 

Xi = — sin(Kr) + Bq . (33) 
r 



with 



Differentiaton gives 

d2C^^B2 
dr r 



= — {Krcos(Kr) — sin(Kr)} . (34) 



3.3 Exterior 

This case differs from the previous one by the substitution /t ^q- Analogously 
we receive the solution 

Xe = — - cos(Kr) H sin(Kr) + Aq (35) 

r r 

with 

Aq = — ^ (36) 

{Ai and A2 constants of integration). For physical reasons this formally correct 
result has to be changed, because it is not in accordance with the boundary condition 
(pTIp): We know that the global source /ic has already be taken into account by 
solving the cosmological differential equations. Therefore in this local context we 
have to drop the inhomogenious contribution by setting Aq — 0, i.e. 

Xe — — cos(Kr) + — sin(Kr) . (37) 
r r 

Differentiation leads to 

^=(A2n-^] -(a,.+ ^) . (38) 

dr \ r J r \ r J r 

In performing these calculations we realized that following basic functions play an 
important role: 

a) [/(z) = z sin z + cos z and b) y(z) = sin z — z cos z . (39) 
Figures |l| and || show the plots of the functions U (z) and V{z). 
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UUzL 




Figure 1: Dark matter basic function U 



VUzL 




Figure 2: Dark matter basic function V 



3.4 Continuity conditions and the final form of the poten- 
tials and of the dm-particle number densities 

At the surface of the spherical central body following continuity conditions have to 
be fulfilled: 

a Xt(ro) ^ Xe{ro) , b — ; . (40 

dr dr 

If these conditions are satisfied, then because of (|20| ) the analogous conditions for 
the dm-particle number densities result: 

a) n.M = n.M, b) ^^^^ (41) 

dr dr 

i.e. thus the continuity of the dm-particle number density is guaranteed, too. Look- 
ing at (|3^) and ( p7| ) we recognize that we have to cope with the 3 constants of 
integration B2 , Ai and A2. From both the conditions (40) we obtain the relations 

a) Ai^-^V{Kro), b) B2 - = - — U{Kr„) . (42) 

K K 

For physical reasons it seems to be sensible to accept that the potential in the whole 
space (interior and exterior) may be determined solely by the central mass. This 
means setting A2 — 0. Then finally the decisive potentials ( ^3|) and (37) including 
their derivatives (E3) and (Isl) read: 



a) Xi = Boil - —U {nro) sin(Kr) 
Kr 



dr Kr^ 
a) Xe = --V{Kro) cos(Kr) 



(43) 



(44) 



dXe _ Bo_ 
dr Kr^ 

Let us mention that in the center of the central body the equations 

a) x^ir = 0)=Bo{l-U{Kro)}, b) ^2^1^-^ = Q (45) 
are valid. 

Furthermore, according to ( |2^ ) we also present with the help of (^3|a) and (^Ja) 
the final form of the distribution functions n,; and Up: 



nmBo f 1 / \ / \ 
a) = 1 - — C/(Kro)sin(Kr) 



Kr 



(46) 



b) ne = ^^^° V{Kro) cos(Kr) . 
ki Kr 



3.5 Series expansions 



For various applications it is convenient to know the first terms of the series expan- 
sion of some functions presented above: 



a) U{z) = l + lz'-lz\ 

b) Viz)^^z'il~^z'): 



(47) 
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Figure 3: Dependence of the interior potential on radius R 



1 2 



b) 



b) 



1 / r 

3 



Vof l + -(-)^--^- 
4^ 1 3Vo^ 15 Uo 



dr 



2 I 5 To 



a) Xe = SoK^r^ 



dXe ^ J_ 

d?" 3r' 



(49) 



(48) 



X*('^ = 0) = --So(Atro)2 



1 - ji^rof 



(50) 



One should realize that the special case k = of ([48|a) and ([49|a) corresponds to 
usual Newtonian physics with the results 



a) XNi = 

b) XNe = 



! " 5 








2ro 





(51) 



1 

3^ 



Concluding this section, we point to the Figures |^ and ^ which show plots of the 
radial dependence of the potential % on r in different ranges (qualitative presentation 
for Bq — 1). Analogously the Figures ^ and || exhibit the radial dependence of the 
dm-particle number density n on r (again qualitative presentation for Bq = 1). For 
these Figures the value k = 0.04[(kpc)~^] was chosen which apparently seems to be 
appropriate for the rotation curves of the Galaxy. 



4 Pressure in the central body 

For physical reasons the radial course of the pressure is interesting, too. Approach- 
ing this aim, here instead starting from the equation of motion of a test body (^) , we 
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Figure 4: Dependence of the exterior potential on radius R 
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Figure 5: Dependence of the interior dm-particle number density on radius R 
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Figure 6: Dependence of the exterior dm-particle number density on radius R 



have to start from our generalization of the Navier-Stokes equation of a continuum 



I dv , c , dhi(Tc\ / ^ V X B . 

h grad <£> H grad s + v = p \ E -\ — gradp — 

dt (Jc dt I \ c ' 



(52) 



-77 ( Ad+- graddivi; 



{p pressure, p electric charge density, rj viscosity). Specialization to the case v = 
(hydrostatics), p = (electro-neutrality), 7] = (perfect fluid) and mass homogene- 
ity leads to the equation 



grad I $ + — s +—1=0 



(53) 



a) grad x H = , 

Mo, 



{p constant of integration). Using ([43|a) we find 



i.e. b) p=-pqx+P 



p = pqBq < — U{K,ro) sin(Kr) - 1 + p . 



(54) 



(55) 



The postulate of vanishing pressure at the surface of the sphere, p{r = vq) = 0, 
leads to the result 



P = MoSoC^(Kro) 



sin(Kr) sin(Kro) 



uro 



(56) 



with 



p(r = 0) = p,oBoU{Kro) 



sin(Kro) 



Specialization to the Newtonian case k = gives the result 



PN 



2ro 



(57) 



(58) 



5 Radial drift of a test body ("Pioneer effect") 

Recent evaluation of radio metric data from Pioneer 10/11, Galileo and Ulysses 
spacecraft indicates an anomalous, constant acceleration acting on the spacecraft, 
being about 50 AU far from the Sun, with a magnitude « 8.5 • 10~^cm/s^, directed 
towards the Sun (Anderson et al. 1998). Apparantly the explanation of this new 
effect seems to be related to a correction of Newtonian mechanics by additional 
terms. But also general-relativistic explanation approaches are under discussion 
(Rosales and Sanchez-Gomez 2000). In the following we try to find an explanation 
of this radial drift effect by applying our theory. 

The equation of motion of a test body (||) with the help of (^a) takes the form 

^^""^ " (59) 



— + gradx + 'i^- , , 
at dt 



0. 



From this equation we learn that for a spherically symmetric potential, as investi- 
gated above, the appearing acceleration A of the body considered can be splitted 
into the Newtonian part and the scalaric part As according to 



A = ^ = A, 
dr 



As. 



(60) 
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Figure 7: Dependence of the additional radial acceleration on radius R 



Referring to the interior and exterior of the central body we identify 



a) A 



Ni 



and 



b) A 



Ne 



(61) 



(internal and external Newtonian accelerations). With the help of (^^), ^8|b) and 
(44d), (^9|b) the additional scalaric parts take the shape 



a) Ast 



b) Ase 



3U{Kro)V{Kr) 



W{Kro)U{Kr) 



1 



-IN Metier 



1 - 



(62) 



'-{nrf 



Remembering that this series expansion is valid for (^r)^ <C 1, we apply this latter 
formula to the empirical facts mentioned above: 



1 - -(nr)^ 



= 8.5 • lO^^cnis 



(63) 



Inserting the well-known values of the Newtonian gravitational constant = 
6.68 • lO^^g^^ cm'^s^^, the mean distance of the spacecraft considered from the Sun 
r = 50 AU, and the mass of the Sun Mc — 1.989 • 10'^^ g, we arrive at the value of 
K for the Sun (1 AU = 1.496 • lO^^ cm) 



K = 3.58 • 10"^^cm"i = 5.36 • 10"''(AU)" 



(64) 



Though our theory predicts a rather constant value of the radial acceleration far 
away from the Sun, nevertheless this acceleration value decreases at even greater 
distance. Figj^ shows the course of the curve in detail. 



6 Dark matter induced new interaction and the 
"5th force" interaction 

Nearly twenty years ago intensive discussions on the validity of the equivalence 
principle and on an eventually existing new interaction, acting also in macrophysics, 
took place. Many proposals were published, where mainly the concrete approach of 
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Fischbach et al. met considerable attention (Fischbach et a. 1986) who introduced 
the notion "fifth force" with the Yukawa type potential energy between two particles 
TOi and 7712 : 



V{r) 



1 + a cxp ( - ^ 



(65) 



(Goo gravitational constant at large distance, a and A free parameters). 

According to our external potential (|9|) the potential energy between such two 
masses considered reads {'^n Gc 



= 7 -^V^(Kro)cos(Kr) « 



Since the functional structures of the potential energies (|6|) and ( |66| ) are dif- 
ferent, the free parameters a , A and k, tq don't correspond uniquely. Roughly 

speaking, only the identification — — is possible. 

A^ 
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